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1 Introduction 

ATLAS [1] and CMS [2] have recently presented the combination of all their Standard 
Model Higgs searches with the full 2011 dataset. Both experiments have registered a 
tantalizing excess of events at about 125 GeV with respect to the hypothesis of the absence 
of the Higgs. However the evidence is not yet sufficient to claim that the Higgs boson has 
been detected and statistical fiuctuations might well explain the data. Therefore, the 
mechanism of Electro Weak Symmetry Breaking (EWSB) remains unclear and the issue of 
high energy vector boson scattering continues to play a central role, either as the final test 
of the nature of the Higgs boson or, if the Higgs doesn't show up, as the main hunting 
ground for clues to alternative explanations. The vector vector scattering amplitudes grow 
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with energy when the bosons are longitudinaUy polarized and violate perturbative unitarity 
at about one TeV [3-8] , requiring either the Higgs or some new physics in the energy range 
accessible to the LHC in order to tame this unphysical behaviour^. 

Many alternative mechanisms of EWSB have been explored. We will not try to sum- 
marize the different models and simply refer to the literature. These theories typically 
predict the presence of new states which, much like the Higgs boson does, keep the scat- 
tering amplitudes small and the full theory amenable to a perturbative treatment. These 
additional particles, if light enough, could be observed at the LHC. If no such state is 
present, the scattering amplitudes become strong as the energy increases and perturba- 
tion theory breaks down. Moving beyond the perturbative approach it is possible to use 
Effective Field Theory (EFT) methods, in particular the Electroweak Chiral Lagrangian 
(EWChL) [13-19], to describe the low energy behaviour of the scattering amplitudes. The 
EWChL is a powerful approach for treating the low energy dynamics of systems with bro- 
ken symmetries. It provides a systematic expansion of the full unknown Lagrangian in 
terms of the fields which are relevant at scales much lower than the symmetry breaking 
scale and does not require a detailed knowledge of the full theory. It is then possible to 
apply Unitarization Methods, using the lowest order terms in the scattering amplitudes 
as building blocks of all order expressions which respect unitarity and agree up to a finite 
order with the perturbative result. 

Scattering processes among vector bosons have been scrutinized since a long time 
[20-34] with an increasing degree of sophistication. In a series of papers [35-38] we have 
analyzed all vector vector scattering channels, 'ijiu, Ajt^f~ , 2j3iu, 2j£^l~vh', 2j£^i^vu 
and 2j4:i (l = e, //) which are observable at the LHC, including all processes at order 
0{a%j^j), 0{a'^j^,ja1) and 0{a'^j^ja'^) as well as tt+jets when appropriate. We have system- 
atically compared a typical SM light Higgs scenario with the Higgsless case and with an 
example of Strongly Interacting Light Higgs (SILH) model [18]. 

In the last few years QCD corrections to boson-boson production via vector boson 
fusion [39-42] at the LHC have been computed and turn out to be below 10%. Recently, 
VBFNLO [43], a Monte Carlo program for vector boson fusion, double and triple vector 
boson production at NLO QCD accuracy, limited to the leptonic decays of vector bosons, 
has been released. 

The first results for the NLO corrections to TV + 4j production have started to appear 
[44] . New techniques which exploit the angular distribution of vector boson decay products 
to determine the ratio of longitudinal and transverse polarization have been proposed in 
[45]. 

A number of papers [46-53] have explored the possibility that the increasing strength 
of the interaction between vector bosons as the energy of the scattering is augmented might 
lead to the formation of resonant states in analogy to Low Energy QCD in which the growth 
of the pion-pion scattering amplitude is regulated by the appearance of the p resonance. 
Most of these efforts have resorted to the EVBA, since strictly speaking the unitarization 
procedure is defined only for on-shell scattering amplitudes between longitudinally polar- 

^ Detailed reviews and extensive bibliographies can be found in Refs. [9-12] 
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ized vector bosons. This approach however suffers of all the well known deficiencies of the 
EVBA and in particular neglects the contribution of transversely polarized vector bosons 
and all off-shell effects which can be sizable. The strong gauge cancellations between 
signal and irreducible background and the reliability of the Equivalent Vector Boson Ap- 
proximation (EVBA) have been studied in [54]. A preliminary analysis in the EVBA of the 
observability of partial unitarization of longitudinal vector boson scattering in SILH mod- 
els at the LHC can be found in Ref. [55]. Unitarization models can describe the strongly 
interacting symmetry breaking sector, with or without formation of resonances. A recipe 
for embedding the unitarized amplitudes in the full, off-shell, 2 — )• 6 amplitude, which is 
necessary for a reliable description of VV scattering at the LHC, has been first proposed 
in Refs. [56, 57] and later refined in Ref. [58] and made available in WHIZARD [59]. 

In this paper we present the implementation of several unitarization schemes, the K- 
matrix, the Inverse Amplitude Method (lAM) and the N/D procedure, within the PHANTOM 
[60] Monte Carlo event generator along the lines of Ref. [58]. The unitarized amplitudes are 
extended off-shell in a natural way and embedded in the framework of a complete six parton 
final state calculation. We then study the prospect of detecting signals of these unitarized 
models at the LHC taking into account all relevant backgrounds, including top-antitop 
production, possibly with additional jets, and vector boson plus four jets production. 

The paper is organized as follows. In Sect. 2 we present the basic formalism of the 
EFT approach and we introduce the unitarization schemes which are available in PHANTOM. 
In Sect. 3 we compare the results obtained in a complete calculation with those obtained 
in the EVBA. Then, in Sect. 4 we present the main results for all relevant decay modes of 
the final state bosons and finally we state our conclusions. 

2 Implementation of unitarized models in PHANTOM 

Pion-pion scattering has been for a long time described by effective Lagrangians. In this 
framework, when the energy increases, non perturbative effects of QCD dynamics are 
unavoidable as for example manifested in the appearance of new resonances. A purely 
phenomenological approach to describe this physics is given by Unitarization Models. They 
are based on ad-hoc formulas that force the amplitudes of Goldstone bosons scattering to 
satisfy the unitarity condition and maintain the low energy behavior. Unitarization models 
are intended to represent the approximate magnitude of these amplitudes beyond the low- 
energy regime and have been able in some cases to describe some resonances of QCD. 

In view of the great similarities between low energy QCD and Electroweak physics, the 
ideas of unitarization models have been translated to a strong symmetry breaking sector 
in several studies [31, 46, 48, 50, 56-58]. They are not complete quantum field theories 
and in particular they typically violate crossing symmetry, but despite these deficiencies, 
unitarization models fulfill their phenomenological purpose of estimating the magnitude 
of strong VV scattering cross sections much beyond the range of validity of the effective 
theory. 

In this section we describe the main aspects of the implementation of unitarization 
models in a complete six-fermion in the final state framework for the PHANTOM Monte Carlo 
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generator. 



2.1 The Electroweak Chiral Lagrangian And Low Energy Amplitudes 

If EWSB is driven by new strong dynamics at TeV scale, the EWChL (inspired by the Chiral 
Lagrangian of QCD) provides the most economical description of electroweak physics below 
this scale. The EWChL accounts for all particles of the SM apart from the Higgs boson. 
The gauge symmetry SU{2)l x U{1)y is maintained by explicitly introducing the three 
Goldstone bosons, uj"'{x) with a = 1,2,3, gathered in an SU{2) matrix field 

E(x) = exp ^ , (2.1) 



where ex" arc the Pauli matrices and v ~ 246 GeV is the decay constant of the Goldstone 
boson that gives the right masses to vector bosons. S transforms under SU{2)i x U{1)y 
as 

S ^ ULix)ml{x), (2.2) 

Vl{x) = exp ), Uy{x) = exp { ^"^y ). (2.3) 

The familiar pieces of the chiral Lagrangian, that emerge for example from the Mh — > 
GO limit of the SM, are: 

£ ='^'IV[(D^E)t(i?''S)] - ^G^.G'^-''^ - ^Wl^W^^^' - \b^.B^'^ 
+ i QlPQl + i QrPQr + i LlPLl + i LrPLr 

- [Ql^MqQr + Ll^MlLr + h.c). (2.4) 

The first term has the form of a non-linear sigma model, which is non-renormalizable in 
four dimensions; therefore, a cut-off scale A must be set. Deviations from the SM in the 
absence of the Higgs boson can be parametrized in terms of a low energy expansion in 

E/ A., consisting of operators of increasing dimension. 

Besides the usual spontaneous breaking pattern of the SM, SU (2) l x U{1)y — )■ U{l)emi 
experiment demands that the Higgs sector also approximately respects a larger, SU{2)i x 
SU{2)r symmetry, S — ^ Ui'^uj^, which is spontaneously broken to the diagonal subgroup 
SU{2)d (where Ul = Ur) by the VEV (S) = 1. Only two dimension-4 operators which 
respect these symmetries are relevant for the study of Vector Boson Scattering (VBS), they 
are: 

A = a4'IV[y^,F^]2, (2.5) 
C^ = a^T^[V^,Vi'\\ (2.6) 

where = (Z)^S)St. 

For the description of the elastic scattering of longitudinal bosons, we are going to 
make use of the Goldstone Boson (GB) amplitudes, which can be translated into the corre- 
sponding physical longitudinal boson scattering through the Goldstone Boson Equivalence 
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Theorem (GBET) [6, 61-63]. The GBET states that an amplitude involving longitudinal 
vector bosons, Vl, is well approximated at high energy by the corresponding amplitude 
obtained by replacing Vl with the Goldstone bosons, uj, in any renormalizable gauge. 
Assuming isospin custodial symmetry, all 2 — )■ 2 Goldstone boson scattering processes can 
be described by a single master amplitude, A{s,t,u) which satisfies A{s,t,u) = A{s,u,t). 
It can be identified as the amplitude of the uj'^ljj~ — t- zz process. 

The lowest order term of the master amplitude in the E/A expansion, identified as the 
Low Energy Theorem (LET) [64] and reproduced by the non- linear sigma model term in 
the Lagrangian (first term in the r.h.s. of Eq.(2.4)), is given by 

AW(s,t,u) = A^. (2.7) 



V 



At next-to-leading order in the E/A expansion, we must include the higher-dimension 
operators, Eqs.(2.5-2.6), and one-loop diagrams. The NLO amplitude is given by [65] 



-(-91/ X , t'^ + u'^ 1 
A^^>{s,t,u) = 4a4 3— + 8a5^ + 



H In 



t{s + 2t) u{s + 2u) f ^j? 



+ 7—^ In — + ^ — - In 

Qv^ \-t 6t;4 



-u 



(2i 



where n is the renormalization scale. The infinities that appear when computing one loop 
diagrams are absorbed by defining renormalized parameters in the higher-dimension oper- 
ators of the effective Lagrangian. The terms proportional to factors of the weak coupling, 
g, which do not grow asymptotically with s, will be recovered in the complete 2 — >• 6 imple- 
mentation. For the moment, we just need the leading contribution from the longitudinal 
boson scattering. 

The 04 and parameters contribute to the T-parameter [66-68] and therefore are 
constrained by electroweak precision data [69, 70]. Stronger limits are obtained if quadrat- 
ically divergent terms are taken into account [71]. Arguments based on unitarity and 
causality could also constrain these parameters [72]. In this case, the magnitude of 04 and 
as is required by data to be smaller than about 10~^. 

2.2 Unitarization Of Low Energy Amplitudes 

In the EWChL framework, when E approaches A, the perturbative expansion starts to 
loose its predictive power. Moreover, low energy amplitudes of longitudinal vector boson 
scattering violate unitarity at much lower scales £^ ~ 1.2 TeV. To describe the magnitude 
of the cross section much beyond the unitarity violation scale we can unitarize the low 
energy amplitudes. 

In order to perform the unitarization procedure, it is convenient to expand the master 
amplitudes into isospin, /, eigenamplitudes according to ^ 



^These expressions as well as those in Eq.(2.10) can be derived using the Clebsh-Gordan coefficients 
for the coupling of two 7 = 1 representation, exploiting the scalar nature of the interaction Hamiltonian. 
The pion states are defined as: |7r+> = l/^/2(|7^l> +i|7r2>) = -|1,1>, |7r°> =|7r3> = |1,0>, \tv~ > = 

1/V2{\-K1> -i\7V2>) = |1,-1> 
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Aq{s, t, u) = 3A{s, t, u) + A{t, s, u) + A(u, s, t), 
Ai{s, t, u) = A{t, s, u) — yl(n, s, t), 
A2{s, t, u) = A(t, s, u) + A{u, s, t). 



(2.9a) 
(2.9b) 
(2.9c) 



The individual scattering amplitudes can be expressed in terms of the isospin eigenam- 
plitudes as follows: 



A{uj~^u> —7- zz) = -Aq{s, 
3 


t,u) - 


1 . / 

-A2{s,t,u) , 


(2.10a) 


A{uj~^z w+z) = ^^i(s, 


t,u) + 


-A2{s,t,u) , 


(2.10b) 


A{U!^UJ~ — )• UJ~^U}~) = -Aq{s, 
o 


t,u) + 


l-Ai{s,t,u) + lA2{s,t,u) , 
2 b 


(2.10c) 


A(a;+w+ tj+a;+) = A2{s,t 


u) , 




(2.10d) 


A{zz zz) = ^Ao{s, 


t,u) + 


2 

-A2{s,t,u) . 


(2.10e) 



Each isospin eigenamplitude is then expanded into partial waves of definite angular 
momenta, J, according to 



Auis) 



d cos 9 Pj{cos6)Ai{s, t, u) 



(2.11) 



where Pj{x) are the Legendre polynomials. 

Notice for instance the behavior of the / = 0, J = eigenamplitude: 



Aoo{s) 



^ {^ai{^l) + lla5(//)) + (21n ^ 



sj 9 V s 



+ 



11 

54 



(2.12) 

The absorptive part of this amplitude is extracted by setting ln(— s) = Ins — ivr, when s > 
and it reproduces the perturbative unitarity relation 



^^-4Xs) = —\Afj{s) 



(2.13) 



where A^^\s) is the lowest order amplitude (the first term of the r.h.s. of Eq.(2.12)) and 
A(2)(s) is the NLO part given by the term proportional to /v^. This amplitude presents, 
up to one loop, the correct singularity structure, crossing property and chiral symmetry. 
Nonetheless, this is not sufficient to guarantee a well behaved cross section at high energies 
such as those of the LHC. 

Therefore, each individual isospin-spin eigenamplitude must be unitarized over the 
whole spectrum employing a specific protocol. We have implemented the three most pop- 
ular models of unitarization: the K-matrix scheme (KM), the Inverse Amplitude Method 
(JAM) and the N/D protocol. The main purpose of these schemes is to transform the 
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isospin-spin eigenamplitudes, Ajj{s), into new expressions, Aij{s), which simultaneously 
respect the unitarity condition, 




and have the appropriate low energy behavior. 



Auis) 



> Aij{s) . 



(2.15) 



2.2.1 K-Matrix Scheme 



The unitarized amphtude through the K-Matrix scheme is given by 




327r 



(2.16) 



327rRe(l/^7j(s)) 



— I 



This expression exactly satisfies elastic unitarity. As an example, let us take the pure Low 
Energy Theorem expression, Aqo{s) = 2s /v'^. The KM procedure transform this expression 
into 



which instead of growing quadratically with energy, asymptotically approaches unitarity 
saturation. 

We have also implemented in the code the possibility of incorporating the exchange 
of heavy resonances of definite angular momentum, J, and isospin, /, unitarized by the 
K-Matrix scheme. We have followed the prescription given in [58], adding new degrees of 
freedom to incorporate resonances coupled to the Goldstone bosons. This unitarization 
procedure normally reproduces resonances with s-dependent widths. 

2.2.2 N/D Protocol 

Besides violating crossing symmetry, the KM unitarization procedure spoils the singularity 
structure of the fixed-order amplitudes. It fixes the absorptive part for s > guaranteeing 
unitarity at arbitrary energy at the cost of ruining the left-hand cut. In the N/D protocol, 
unitarity is exactly restored with the extra quality of improved analytical properties. In the 
N/D method, each partial wave amplitude is expressed as the quotient of two functions. 



The denominator function, D{s), contains the right hand cut (or unitarity cut) while the 
left-hand cut is incorporated in the numerator function, N{s). 

In Eq.(2.12) for instance the left hand cut first appears through the term (7/9) ln(/i^/s) 
which acquires an imaginary part for s < 0. To unitarize the amplitudes and simultaneously 
reproduce the left hand cut of Eq.(2.12), at one loop precision, we follow the prescription 
given by [49] . We define 




(2.17) 




NiAs) 



(2.18) 




(2.19) 
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where M is a new free parameter. The unitarized amphtudes take the form 

Nrj{s) = 4]{s) + Af]{s) + G{s){4]{s))' ; (2.20) 
Dij{s) = l + Gis)Njjis). (2.21) 

The terms ^jj(s) are originated from the lowest order ampUtude, Eq.(2.7), while the 
AjJ{s) are originated from the NLO terms, Eq.(2.8). Njj{s) bears exclusively the left hand 
cut, which is kept at the unitarized amplitude, A^j^{s), at the one- loop level precision. 



For s > 0, Njj{s) is real and the unitarity cut is completely provided by the denominator 
function, j 
Eq.(2.14). 



function, Djj{s), making Afj^{s) = ^Yj\^s) *° exactly respect the unitarity condition. 



2.2.3 Inverse Amplitude Method 

In the Inverse Amplitude Method (lAM), isospin-spin eigenamplitudes are unitarized by 
the following prescription: 

l-Af]{s)/AY]is) 

It can be shown that Eq.(2.22) is a special case of the N/D method, it maintains the 
proper analytical structure with the correct branch cuts. In addition, for certain values 
of the chiral coefficients, 04 and as, the unitarized amplitudes, both by N/D and lAM 
protocols, present poles that can be interpreted as dynamically generated resonances. 

This method has been widely and succefully applied in pion-pion and pion-kaon scat- 
tering, in many cases reproducing lightest resonances. 

2.3 Implementing Unitarized Models in Six Parton Final States within PHANTOM 

No detailed phenomenological study has been performed in the context of unitarization 
models in a complete calculation. In most cases, previous results have relied on the Effec- 
tive Vector Boson Approximation (EVBA), which is known to produce inaccurate predic- 
tions for VBS at high energy. With a complete 2 — )• 6 calculation, all diagrams are summed 
coherently, interference terms and off-shell effects are completely accounted for. The con- 
tribution of transversely polarized bosons is also correctly included. A detailed description 
of the implementation of K Matrix unitarized model in WHIZARD and some examples of 
simulation with the full 6 fermion calculation is given in Ref. [59]. 

It is important to notice that the LET part is already present in the complete tree-level 
calculation. In order to avoid double counting it is necessary to define a correction to the 
on-shell amplitudes, given by 

AAjj{s) = Aij{s)-A^^]{s). (2.23) 

These corrections to the isospin-spin eigenamplitudes are translated back into isospin eige- 
namplitudes according to 

AAi{s, t) = ^(2 J + 1) AAi J {s)Pj {cos 6) . (2.24) 
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For our purposes, it is enough to truncate the series at J = 2. Finally, the isospin eigenam- 
plitudes are translated into corrections to the individual scattering amplitudes according 
to Eq.(2.10). 

We have for example: 



zz 



+ 4 



24s2 
5v^ 



(A^oo(s) - AA2o(s)) 



5v^ 
12s' 



(AAo2(s) - AA22{S)) 



4s2 



(AAo2(s) - AA22{s)) 



"1 t'^ + 



(2.25) 



In order to introduce each of these on-shell elastic scattering amplitudes into the com- 
plete 2 — )• 6 matrix elements, a method based on the definition of new quartic vertexes 
has been suggested in [56, 57] and further developed in [58]. Our implementation of this 
approach in the PHANTOM Monte Carlo generator is briefly explained in the following. 

2.3.1 New Quartic Vertexes 

In unitary gauge the quartic gauge interactions of massive vector bosons can be extracted 
from the EWChL, Eqs.(2.4-2.6), with S = 1: 

Cqgc = cos^ e„ [gf^Z^'Z''W-W+ - gl^Z^^Z^W-'W^] 



+ 



+ 



9 



2 



5f ^ {W'^W^ 



h^^{Z>'Z^ 



(2.26) 



Considering only the £4 and £5 extra terms, Eqs. (2.5-2.6), the values of the couplings are 
given by [59] 



1 + 



9 



cos4 (9, 



-04 



WW 1 , 2 

9i =1 + 9 a4, 
^zz _ „2 



92 

„WW 
92 



1 



9 



-as , 



cos4 (9^ 
1 - 5^(04 + 205) 



g (04 + as) 



(2.27) 



In the SM, g^^ = 1 and h^^ = 0, where i = 1, 2 and VV = WW, ZZ. 

As already mentioned, the unitarization procedures used herein violate crossing sym- 
metry, therefore, it is not possible to write a Lagrangian for the model and to derive the 
corresponding vertexes in a general form. Instead, it is necessary to identify first which kind 
of vector vector scattering processes are embedded in the complete 2 — )• 6 set of Feynman 
diagrams. As a second step the vertexes are appropriately modified. The possible scatter- 
ing reactions, W+VF" ^ W^W' , W+W' ZZ, W^W"^ W^W^, can be identified 
from the flow of fermion external momenta. 
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For each scattering type, the quartic vertex is modified in such a way to reproduce the 
full on-shell scattering amplitudes. They can be written in the form: 



■ 2 2 

-ig cos I 



ig'^[2kY^g^'Pg'"'- kf^g^'^g'^P- kf^gP^gP"] 



ViZPZ^'ZPZ'') = ig"' 



cos4 ft 



.[^p^M-^P- + ki^gP-g-P + ki^gPPg' 



(2.28a) 
(2.28b) 

(2.28c) 



where the ki coefficients are form factors which depend on the particular vector vector 
scattering process. For example for ZZ — )• W~^W~ we have: 



zw 



1 



cos^ 6^ 



^ (AAoois) - AA2o{s)) - ^ (AAo2(s) - A^22(s)) 



K,2 



..ZW 



1 + 



COS' 



.4 0. 



5v^ 
4? 



(A^02(S) - A^22(S)) 



(2.29a) 
(2.29b) 



Notice that for on-shell VBS, these new vertexes reproduce the unitarized amplitudes (e.g. 



Eq.(2.25)) at high energies. Since 



in the high energy limit the Mandel- 



stam variables of the scattering can be reexpressed in terms of the vector longitudinal 
polarizations in the form: 



s2 = 2pi ■ p2 2p3 ■ p4 ~ AMv,Mv,Mv,Mv, g^,u9pa ^1^24^4, (2-30a) 
= 2pi ■ p3 2p2 ■ PA w ^Mv^Mv^Mv^Mv^ g^,pgua e^2e44, (2.30b) 
= 2pi ■ Pi 2p2 ■ P3 - 4.Mv,Mv,Mv,Mv, g^^aQup 4^24^1^ (2.30c) 

where the subindexes 1, 2 indicate the two incoming bosons and 3, 4 the two outgoing ones. 
My are the masses of the bosons, either My/ or Mz- With this identification, the corrected 
amplitudes take exactly the form of quartic vertexes contracted with external polarization 
vectors. The difference with respect to the complete six-fermions final state calculation is 
that the role of the polarization vectors in Eq.(2.30) is played by the sub-diagrams with 
the final state fermions and the boson propagators. 




(2.31) 



At low energies, the corrections Avl/j tend to vanish; therefore, the complete cal- 
culation is recovered. For high energies, the typical growth of longitudinal vector boson 
scattering is moderated by the A^/j corrections reproducing the unitarized elastic scatter- 
ing of the GB, which is a good approximation to the scattering of the longitudinal modes 
at high energy according to the Goldstone Boson Equivalence Theorem. 
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3 Comparison with EVBA results 



The results of our implementation have been compared with those obtained in Ref. [50] 
by Butterworth, Cox and Forshaw, where an on-shell calculation based on the Effective 
W Approximation has been employed. Three instances of the Inverse Amplitude Method, 
Eq.(2.22), have been considered. They differ in the choice of the 04 and 05 coefficients as 
reported in Tab. 1. The renormalization scale has been set at TeV. 



scenario 


a4(lTeV) 


«5(lTeV) 


A 


0.0 


0.003 


B 


0.002 


-0.003 


D 


0.008 


0.0 



Table 1. Parameters for the lAM models. 



The lAM procedure gives rise to poles in the unitarized amplitudes which can be 
interpreted as resonances whose mass and width are given by Eq.(3.1) for the vector channel 
1 = 1, J = 1 and by Eq.(3.2) for the scalar case 1 = 0, J = 0. Scenario A presents a scalar 
resonance at about 1 TeV, scenario B a vector resonance at about 1.4 TeV, and scenario 
D a scalar resonance at 0.8 TeV and a vector one at 1.4 TeV. 



9(47r)2 

Tor 



Generation cuts 

Pt{^^) > 20 GeV 

\r]{£^)\ < 3.0 
PtU) > 30 GeV 

\v{j)\ < 6.5 

M{jj) > 60 GeV 
M > 20 GeV 

Table 2. Standard acceptance cuts applied in the event generation and present in all results of 
this section. Here j = d, u, s, c, b, g. 



We have computed the complete set of purely electroweak processes for the 2jii e^vv 
final state, with the generation cuts in Tab. 2. A big contribution to these final states 
comes from electroweak ti production. Therefore, in order to avoid this contamination in 
the comparison, we have considered only processes without bottom quarks in the initial or 
final state. Moreover, different cuts in the pseudo-rapidity difference between tag jets have 
been tested in order to enhance the scattering contribution. 
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Figure 1. Comparison between the M{WW) distribution in PP — 2i^,^e^vv in the complete 
calculation and the EVBA calculation in Ref. [50] (in purple) multiplied by the appropriate branch- 
ing ratio, BR = 2/81, for the lAM A model. The complete calculation is shown with different cuts 
in pseudo-rapidity. Processes with external b-quarks have been discarded. 

The VV invariant mass distribution in the lAM A model is presented in Fig. 1. It 
shows the prediction of the complete calculation, as implemented in PHANTOM, and the 
corresponding cross section derived in Ref. [50] in the EVBA. The data points for the EVBA 
results have been extracted from the figures of Ref. [50] and rescaled by the branching 
ratio for W~^W~ — )• ^i^e^vv, BR = 2/81, since the W decays were not included. For the 
complete calculation, four different cuts on the difference in pseudorapidity between jets, 
^vijj) > 3, 4, 5, have been applied. As expected the resonance peak is more clearly 
seen at larger separation between tag jets. 

It can be observed that the curve with Arj > 4 is well reproduced by the approximation 
in Ref. [50] close to the peak at about one TeV. Outside this regions the curve of Ref. [50] 
underestimate by a large amount the actual cross section. The comparison between the 
present results and those obtained in EVBA should not be taken too literally. In EVBA 
the final state jets have been completely integrated over and therefore no separation cut 
between tag jets can be applied. While the approximate total cross section obtained is in 
rough agreement with the full calculation in the peak region, too many details would be 
missing in a simulation based on the EVBA to allow for a realistic study of the experimental 
observability of these processes. 

In Fig. 2 we show the comparison for all three scenarios. As before, Ar] > 4 and no 
process with external b-quarks is included. Again, large discrepancies can be observed 
outside the peak region. We also present in black the result for the no-Higgs scenario. It is 
interesting to notice that the unitarized models and the no-Higgs scenario are in reasonable 
agreement outside of the resonance peak in the considered mass range. We will comment 
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Figure 2. Comparison between the M{WW) distribution in PP — 2j^^e^vv for the lAM A [at 
top left), lAM B (ai top right) and lAM D {bottom) models obtained with a complete calculation 
with Ary > 4 and the predictions from Ref. [50] multiplied by the appropriate branching ratio, 
Si? = 2/81. 

further on this point when dicussing the different final states in Sect. 4. 

4 Counting experiments for unitarized models at LHC 

In this section, a number of phenomenological studies at the LHC for some typical choice of 
unitarization schemes and parameters are presented. We limit our analyses to the design 
energy of 14 TeV and we do not attempt to determine the parameter range for which 
resonances in the new scenarios can be experimentally discovered. For some final states 
and for scalar resonances, looking for new resonances is just an extension of Standard Model 
Higgs searches. For different kind of resonances however a completely different strategy 
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might be required. We will therefore examine in which cases there will be a significant 
excess of events at high invariant mass and if it will be possible to detect it at the LHC. 
For this reason we give to these searches the name of counting experiments. All our results 
are at parton level only, and must be considered as a first indication of the effects of 
unitarized models in the context of a complete lowest order calculation which avoids the 
EVBA and takes all relevant irreducible backgrounds into account. 

In addition to the cases already considered in the previous section (lAM A, lAM B 
and lAM D) we will consider the following scenarios: 

- lAM C: it is another Inverse Amplitude Method unitarization model with parameters 

a4 = 0.002 and 05 = —0.001. This model produces a vector resonance at about 1.9 
TeV. It was already discussed in Ref. [50] . 

- N/D A: uses the N/D unitarization method with M = 1 TeV, 04 = and = 0.003. It 

presents a 1 TeV scalar resonance as in the lAM A model but with a softer resonance 
peak. Notice however that in general the N/D method gives models which can be 
quite different from the other methods we have considered. 

- KM-LET and lAM E: these two scenarios contain no resonances at all. In the first 

case the Low Energy Theorem divergent contribution, neglecting all NLO terms, is 
unitarized using the KM method. In the second case the lAM procedure with 04 = 
and as = is employed. They will be compared to the no-Higgs scenario {Mh — )■ 00 
in the Unitary Gauge SM) where no unitarization is performed. 

A detailed description of the comparison among all scenarios will be given for the 
channel 2jW'^W~ — t- 2ji^i^i'V', while for all other channels we will limit ourselves to the 
lAM A, lAM B, lAM D, KM-LET and no-Higgs cases only. 

We have generated about half a million unweighted events for each channel and scenario 
and for the main backgrounds. All 0{a%j^,j) and 0{a%j^fa1) processes have been included 
in our definition of the signal. At generation level the cuts of Tab. 3 have been used. For 
each specific channel we have tried to determine which additional cuts could improve the 
separation between the unitarized model results and those obtained in the SM with a light 
Higgs . The specific cuts used at analysis level for each channel considered in the following 
are reported in Appendix A. All results reported below are obtained after both basics and 
analysis cuts. 

4.1 2jW+W- 2ji+l-uv 

The 2jl'^(l.~vi> final state is one of the most important channels for the study of alternative 
symmetry breaking mechanisms at the LHC [38]. It is one of the channels in which the 
presence of two neutrinos in the final state does not allow the reconstruction of the invariant 
mass of the boson pair produced in association with two jets. This may appear as a limit 
since the invariant mass of the two bosons represents the center of mass energy in the case 
of the underlying boson boson scattering, and it is therefore the most appropriate variable 
to examine in order to discover signs of strong boson boson scattering. However, these 
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Basic Cuts 

Pt{^^) > 20 GeV 

|t?(^^)| < 3.0 

M{e+r) > 20 GeV 

M{l+l-) > 250 GeV {2jW+W-) 

76 GeV < M{i+e-) < 106 GeV (2jZZ) 

PtQ-) > 30 GeV 

< 6.5 

M(ji) > 60 GeV 

Mjjfh) < 70 GeV; Mjjfj,) > 100 GeV 
|A770j)I > 3.0 (2i2£2i.) 

|A??0/jb)| >4.0 {2jAe, Ajlu, AjU) 

\M{jjj) - Mtop\ > 15 GeV {Ajiu, AjU) 
\M{jiurec) - Mtopl > 15 GeV {3iu + 2j, Ajiv) 
70 GeV < M(jJc) < 100 GeV {4.jlv, ijii) 
ARjjj) > 0.3 i4jiu, AjU) 

Table 3. Basic cuts used at generation level. 



final states are not affected by the huge QCD background which is typical of 4 jets plus 
two leptons final states. It has been shown [33, 37] that using invariant mass distribution 
of the two charged lepton can give very good results in discriminating between the light 
Higgs SM case and other scenarios. 

As obvious, both 2jW+W- 2jl+l-v9 and 2jZZ 2jl+l-u9 contribute to this 
final state when the flavour of the two final state leptons are the same. However, cutting 
on the invariant mass may easily discriminate between the two channels. 

We start considering the case of 2jW^W~ — )• 2j£~^i~i'u and we will use this channel 
to present a detailed analysis of the main features of the various unitarized models we 
consider. 

In Fig. 3 one can see the invariant mass distributions of the WW pair in the different 
models. 

These distributions cannot be measured in practice because the invariant mass of two 
leptons and two neutrinos cannot be reconstructed. It is however possible to compute them 
at MC level and we present them here in order to clearly see the behaviour of the various 
models. 

On the left hand side of Fig. 3 the distributions for the three non resonant scenarios are 
reported and compared with the prediction of the SM with a light Higgs. Both unitarized 
distributions and the no- Higgs one are larger than the SM one and the difference between 
strongly and weakly coupled theories can be easily recognized. It is worth mentioning that 
the separation between the SM result and those for the no-Higgs, KM-LET and lAM E 
cases is sharply increased by the additional cuts in Appendix A which enhance the boson 
boson scattering contribution. The two unitarized models are very much similar, even if 
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M(WW) [GeV] M(WW) [GeV] 



Figure 3. WW mass distributions in the 2jW~^W — > 2j£^£ vv channel for different scenarios. 

the unitarization scheme is different and in the lAM E case the NLO contributions are 
included. The no-Higgs scenario is practically indistinguishable from the previous ones up 
to an invariant mass of about 1.5 TeV. Above this threshold the effects of the violation 
of unitarity in the no-Higgs model manifest themselves. Because of this violation the no- 
Higgs model cannot be a consistent theory and can be valid only below an invariant mass 
of the order of a TeV. It has to be noticed however that at the LHC with design energy or 
lower, the vector vector pairs produced in the bulk of the events will have a mass smaller 
than 1.5 TeV and as a consequence the violation of unitarity will be of small practical 
relevance. Therefore the no-Higgs scenario can definitely be considered a good benchmark 
model to analyze the possibility of detecting new physics signals. This is also evident 
from the right hand side plot in Fig. 3, where one can recognise the peaking structure of 
the various unitarized models we have considered. If the distribution could be actually 
measured, the number of events of the various scenarios which develop resonances would 
be much larger than the ones of the left hand side. It is also evident from the right hand 
side plot of this figure that both scalar and vector resonances manifest themselves, as a 
consequence of the fact that the boson boson scattering subdiagrams in this case correspond 
both to W^W- W^W (which prod uce both a vector and a scalar resonance) and to 
ZZ — )• W^W~ (which produce a scalar resonance). 

In Fig. 4 the distributions of the invariant mass of the two charged leptons are pre- 
sented. These distributions can be experimentally reconstructed and can be employed to 
discriminate among the different physical models. As it is obvious in these distributions 
all peaking structures have been smeared out but it is clear from the figures that the new 
physics models and the benchmark no-Higgs case can be clearly distinguished from a light 
Higgs scenario. The difference between the light Higgs case and new physics grows with the 
// invariant mass for the non resonant models. On the other hand the number of events at 
high invariant mass decreases. One has to find a good compromise in selecting a range of 
masses which guarentees a good statistics and a sizeable difference between unitarization 
models and the SM. 
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M(ll) [GeV] M(ll) [GeV] 

Figure 4. Lepton-lepton mass distributions in the 2jW~^W~ — >■ 2j£~^£~iyp channel for different 
scenarios. 

This behaviour is exploited in Tab. 4, where the cross sections are presented as a 
function of the minimun M[ll) invariant mass. In the table the contribution of the tijj 
background which can be misinterpreted for the signal when the h' s go undetected is also 
reported [38]. 



Mcut 


no-H 


KM 


lAMA 


IAMB 


lAMC 


lAMD 


lAME 


N/D 


SM 


tijj 


300 


.337 


.303 


.631 


.400 


.412 


.867 


.355 


.367 


.179 


.173 


400 


.212 


.186 


.413 


.274 


.277 


.547 


.224 


.240 


.100 


.0890 


500 


.139 


.115 


.246 


.190 


.187 


.304 


.142 


.150 


.0577 


.0407 


600 


.0968 


.0724 


.132 


.130 


.126 


.160 


.0897 


.0931 


.0332 


.0215 


700 


.0696 


.0461 


.0658 


.0862 


.0858 


.0898 


.0571 


.0568 


.0217 


.0138 



Table 4. Total cross section (in fb) for the {W^W )£^£ vv + 2j channel with the full set of cuts 
in Tab. 3 and Tab. 18 in function of the minimum l£ invariant mass, M{U) (in GeV). 



One can study the Probability Distribution Functions (PDFs) for the number of events 
which can be measured at a given luminosity. This exercise can be repeated for various 
Mcut in order to find the optimum cut for which the probability of dicriminating between 
the light Higgs case and new physics scenarios is largest. We show only the PDF's for the 
optimized cut. They are presented in Fig. 5 for two reference luminosities L = 50 fb^^ and 
L = 200 fb"^ Here and in the following, the PDF's are computed assuming Poissonian 
statistical fluctuations of the number of events computed by the MC and a theoretical error 
on the number of signal events which we model as a flat distribution of ±30% from the 
actual value [37, 38]. 

One can notice that the curves corresponding to models with no resonances (KM-LET, 
no-Higgs and lAM E) are closest to the light Higgs distribution and are not much different 
from each other. The ones which are most separated are the ones corresponding to models 
with resonances, the separation being larger for lower mass resonances. In the plots the 
vertical line represents the 95% limit of the light Higgs distribution. We therefore compute 
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what we call the PBSM@95%CL (Probability Beyond the SM at 95% Confidence Level) 
for the various scenarios as the probability that a number of events larger than the 95% 
limit occurs. The results, reported in Tab. 5, show quantitatively that the probability of 
excluding the no-Higgs model lies among the other non resonant ones. Models with low 
mass resonances can be easily excluded or verified already with luminosity L = 50 fb~^. On 
the contrary a higher luminosity is required for models which are non resonant or contain 
high mass resonances. 



L{ fb-i) 


Mcut 


no-H 


KM 


lAMA 


IAMB 


lAMC 


lAMD 


lAME 


N/D 


50 


300 


49.4% 


37.8% 


97.2% 


68.5 % 


71.6% 


99.9% 


55.4% 


59.2 % 


200 


400 


82.5% 


68.4 % 


100% 


97.2 % 


97.5% 


100% 


87.2% 


91.8 % 



Table 5. Probability to find a number of events larger than 95% limit of the SM (PBSM@95%CL) 
for the various scenarios in the 2jW+W~ — >■ 2jt^i~vv channel for two different luminosities L = 
50 fb"^ and L = 200 fb~^ M^ut is expressed in GeV. 



4.2 The other two neutrinos final states 

The final state with two jets, two neutrinos and two charged leptons of different sign can 
result from WW or ZZ production and decay. In a complete six fermion approach, when 
the two leptons are of the same flavour, the separation of the two channels has to rely on 
cuts on the mass. Another well known and interesting channel for studying boson 

boson scattering is 2jl^£^vv in which the two leptons are of the same sign. These two 
channels will be discussed in the following in the context of unitarized models. 

4.2.1 The 2jZZ 2ji+i-uu channel 

As already mentioned, this channel has been separated from the 2jWW — >• 2jl'^i~vv with 
same flavour leptons requiring \M{U) — Mz\ < 15 GeV. A ZZ final state can be produced 
in two different scattering processes: ZZ — t- ZZ and W~^W~ — )• ZZ. As a consequence 
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this final state has many analogies with the channel we have previously examined. For 
such a reason we do not present here the final state (not measurable) ZZ invariant mass 
distribution. Instead we consider the transverse ZZ mass: 

Ml{ZZ) = [^MI+pI{U) + ^Af| - \pTm+PTmiss\^. (4.1) 

The transverse mass distributions for the lAM A, lAM B, lAM D and KM-LET models 
after basic and analysis cuts are reported in Fig. 6 and compared with the no-Higgs and 
the light Higgs SM scenarios. It is rather evident that with the full set of cuts the non 
SM scenarios differ significantly from the light Higgs case. It is in particular interesting 
that in this channel the lAM A and lAM D scenarios are the two models which differ the 
most from the SM predictions. They both present a relatively light scalar resonance. On 
the contrary for the lAM B model the signal is not so prominent since its vector resonance 
cannot appear in this channel. 



Mt(ZZ) 




400 600 800 1000 1200 1400 1600 ' 1800 ' 2000 
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Figure 6. Transverse mass distributions of the ZZ-system in the IjZZ — ^ 2j£'^£ vv channel with 
full set of cuts, Tab. 3 and Tab. 18 of Appendix A. 

Even if the transverse mass does not directly correspond to vector vector scattering 
center of mass energy, the KM-LET and no-Higgs models have, as before, a similar be- 
haviour for low invariant masses while they start to deviate from each other above 1 TeV. 
In terms of cross sections their difference is only of the order of 10% if we do not restrict 
ourselves to very high invariant mass cuts. The results for the cross sections are reported 
in Tab. 6 . The expected number of events is rather small. 

From the values reported in the Tab. 6, we have determined the cut which maximizes 
the separation of the different scenarios which turns out to be Mcut = 500 GeV for L = 
50 fb~^ and Mcut = 600 GeV for L = 200 fb~^. The corresponding Probability Distribution 
Functions are shown in Fig. 7 and confirm that for the lAM A and lAM D models we are 
nearly certain to be able to observe a deviation from the SM already with 50 fb^^. For 
all considered models, with 200 fb^^ of data, the probability of observing a discrepancy 
between the experimental results is above 80%. 

If one compares Tab. 7 with Tab. 5, it is apparent that in order to discover signals of 
new physics in boson boson scattering processes, one must examine all possible channels 
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J^cut 


no-Higgs 


KM— Lilj i 


T A i\ /r A 
1AM A 


T A 1\ /r T3 

lAM a 


T A A /r 

lAM U 


dM 


300 


.143 


.125 


.249 


.109 


.319 


.0540 


400 


.120 


.105 


.224 


.0886 


.291 


.0396 


500 


.0887 


.0739 


.181 


.0589 


.233 


.0214 


600 


.0691 


.0537 


.145 


.0408 


.172 


.0118 


700 


.0547 


.0395 


.110 


.0285 


.103 


.00697 



Table 6. Total cross section (in fb) for the {ZZ)t'^t vv + 2j channel with the full set of cuts in 
Tab. 3 and Tab. 18 in function of the minimum Mt{ZZ) (in GeV). 



PDFs PDFs 




Figure 7. PDFs of the number of events in the {ZZ)£+£^iyv + 2j channel for Mcut = 500 GeV, 
L = 50 fb^i (left) and M^ut = 600 GeV, L ^ 200 fb"^ (right). 



because some models can produce a modest signal in some channels while being easily 
detectable in others. 





Mcut (GeV) 


no-Higgs 


KM-LET 


lAM A 


lAM B 


lAM D 


50 


500 


68.8428% 


56.5684% 


97.3776% 


41.4185% 


99.4127% 


200 


600 


98.5618% 


93.7842% 


99.9998% 


80.8113% 


100% 



Table 7. PBSM@95%CL, defined in Sub-Sect. 4.1, for the various scenarios in the 2jZZ 
2jl^l^vv channel. 



4.2.2 The 2j(.^l^vv channel 

The channel with two neutrinos and two same sign leptons in the final state has long been 
considered one of the most promising because of the absence of large backgrounds. We 
will however see that the relevance of this channel depends quite a bit on which theory we 
are examining, as this final state corresponds to an underlying scattering of two same sign 

As usual we start considering the invariant mass distributions of the two leptons. 
They are shown in Fig. 8. The full set of cuts has been imposed. No peaking structure is 
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present. This reflects the fact that, for the models we have examined, neither scalar nor 
vector resonances can be produced in this channel. In this case the no-Higgs and KM-LET 
models, which do not develop resonances in any channel, give the largest predictions. The 
lowest curve is as usual the SM light Higgs scenario. All the others fall in between. 



M(ll) 




200 400 600 800 1000 1200 1400 1600 1800~ 2000 

M(ll) [GeV] 



Figure 8. Distributions of the tl invariant mass for the 2it^l^vv channel with full set of cuts, 
Tab. 3 and Tab. 18. 

These results are made quantitative by the values of the cross sections as a function 
of the dilepton invariant mass cut which is shown in Tab. 8 and by the PDF distributions 
of Fig. 9. 



Mcut 


no-Higgs 


KM-LET 


lAM A 


lAM B 


lAM D 


SM 


200 


.435 


.407 


.310 


.332 


.254 


.206 


300 


.290 


.267 


.183 


.201 


.141 


.114 


400 


.191 


.171 


.107 


.120 


.0768 


.0629 


500 


.129 


.112 


.0643 


.0744 


.0430 


.0351 


600 


.0886 


.0760 


.0403 


.0474 


.0250 


.0194 


700 


.0614 


.0517 


.0250 


.0304 


.0150 


.0112 



Table 8. Total cross section (in fb) for the t^V^vv + 2j channel with the full set of cuts in Tab. 3 
and Tab. 18 in function of the minimum invariant mass, M{ll) (in GeV). 

One sees immediately that the number of events produced in this channel is not large. 
Moreover the lAM A, lAM B and lAM D models can hardly be distinguished from the 
light Higgs SM case. Tab. 9 shows however that non resonant unitarized models can give 
detectable signals of new physics at high luminosity. 



L{ fb-i) 


Mcut (GeV) 


no-Higgs 


KM-LET 


lAM A 


lAM B 


lAM D 


50 


300 


81.42% 


74.36% 


34.48% 


44.36% 


13.60% 


200 


400 


98.96% 


96.92% 


57.72% 


72.03% 


17.21% 



Table 9. PBSM@95%CL, defined in Sub-Sect. 4.1, for the i^t^w + 2j channel. 
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In conclusion, one has to be very careful in relying on this channel for detecting new 
physics because, even if it has the advantage of very small backgrounds, it is not sensitive 
to the types of resonances considered here. On the other hand it achieves one of the best 
dicriminating powers for scenarios without resonances. 

4.3 Final states in which the boson boson mass can be reconstructed 

The final states in which the boson boson mass Myy can be reconstructed are those in 
which at most one neutrino is produced, namely + 4j, + 4j, 3£i^ + 2j and Al + 2j. 
Reconstructing Myy is very useful because it corresponds to CM. energy of the underlying 
boson boson scattering. As it happens for on-shell boson boson scattering, the differences 
between weakly and strongly interacting theories manifest themselves more clearly at high 
invariant mass. 

The Al + 2j channel, with two opposite sign charged lepton pairs has a very low cross 
section and it corresponds to a ZZ final state. 

In order to compute the invariant mass of the two bosons, when a neutrino is present 
in the final state its longitudinal momentum is reconstructed with the usual procedure of 
forcing the invariant mass of the ^v pair to be equal to the W boson nominal mass, 

{p'+p^f = M^^. (4.2) 

This determines, up to a two fold ambiguity, the longitudinal component of the neutrino 
momentum [36]. 

The Aj^l' and + 4j channels contain 4 final state jets. Two of them, the most 
forward and most backward one, will provide tagging and will be required to be well sepa- 
rated and energetic. The two central ones will be considered as candidates to reconstruct 
an electroweak boson. The boson boson invariant mass will be assumed to be that of the 
system formed by the two central jets, the charged lepton and the reconstructed neutrino. 

The cross section of the four jet channels is relatively large, however they are affected 
by substantial backgrounds. In particular V + 4j QCD processes give very large rates 
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even after cut optimization. For the Ajiv final state also tt and tt + jets production 
has to be considered carefully. We have discussed in detail in our previous papers [36- 
38] how to deal with these backgrounds. Here we only recall that when computing the 
Probablity Distribution Functions we will attribute a statistical error to these contributions 
but no theoretical error since we assume that they can be measured in nearby regions and 
extrapolated to the kinematic range of interest for VV scattering. 

4.3.1 The Aj^u channel 

The final state of this cannel corresponds to both ZW and WW production. The invariant 
mass distributions of the boson pair with the full set of cuts is reported in Fig. 10. The 
corresponding cross sections for various values of the invariant mass cut are given in Tab. 10. 
Both scalar and vector resonances are clearly visible in this channel. The W+Aj background 
is much larger than the SM electroweak prediction; it turns out to be of the same order 
as the expected yield of the KM-LET and no-Higgs non resonant models. It must be 
pointed out that this background, computed at tree level with MADEVENT at a fixed scale, 
has been rescaled to the scale Q = -^^/2, H'j. = YliPrU) + Et{W), where E^{W) = 
Pj^{W) + M'^{W) [44]. With this dynamical scale the W + 4j contribution has sensibly 
decreased compared to the fixed scale results. 
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Figure 10. Distributions of the mass of the FW^-system for the i^i' + Aj channel with full set of 
cuts, Tab. 3 and Tab. 18. 



Mcut 


no-Higgs 


KM-LET 


lAM A 


lAM B 


JAM D 


SM 


W + Aj 


tt + 2j 


600 


2.36 


2.23 


3.66 


3.04 


5.46 


1.048 


2.03 


.432 


800 


1.558 


1.46 


2.56 


2.32 


3.36 


.618 


1.15 


.167 


1000 


.966 


.877 


1.13 


1.76 


1.90 


.338 


0.62 


.0617 


1200 


.526 


.478 


.414 


1.26 


1.27 


.188 


0.30 


.0264 



Table 10. Total cross section (in fb) for the 1^1/ + 4j channel with the full set of cuts in Tab. 3 
and Tab. 18 in function of the minimum jcjc^i^ invariant mass (in GeV). 
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Figure 11. PDFs of the number of events in the i^i^ + Aj channel for Mcut = 800 GcV, L = 50 fb"^ 
(left) and Af^t = 800 GeV, L = 200 fb"^ (right). 



Examining the PDF distributions of Fig. 11, one realizes that the situation is sim- 
ilar to the one already discussed for the 2jW~^W~ — )• 2ji~^£~ui' channel: the models 
with the lightest resonances are the ones which can be most easily distinguished from the 
Standard Model with a light Higgs. The KM-LET has the smallest cross section of all uni- 
tarized models. It is however well separated by the SM case especially with a luminosity 
of 200 fb-^ The no-Higgs model behaves similarly to the KM-LET one, confirming that, 
even if unphysical, this model has good properties as a benchmark case. 



L( fb-i) 


Mcut (GeV) 


no-Higgs 


KM-LET 


lAM A 


lAM B 


JAM D 


50 


800 


94.51% 


91.03% 


99.99% 


99.97% 


100% 


200 


800 


99.93% 


99.64% 


100% 


100% 


100% 



Table 11. PBSM@95%CL, defined in Sub-Sect. 4.1, for the + 4j channel. 

Finally, the number of events is rather large already at L = 50 fb^^ and the discrimi- 
nating power described by the PBSM@95%CL values of Tab. 11 is rather encouraging. 

4.3.2 The l+i' + Aj channel 

This channel is in many respects similar to the one we have just considered. The final state 
in this case corresponds to ZW or ZZ production in the boson boson scattering processes. 
There are no contributions from ti + 2j but the Z + 4j background is large. It has been 
evaluated with the running scale Q = Hip/ 2 as before. 

In Fig. 12 we show the distribution of the invariant mass of the two final state bosons 
for the different models. After all cuts, the signal cross sections in this channel are about 
five times smaller than those of the 4:j£v channel, as shown in Tab. 12. This implies that 
the discriminating power for new physics models is sensibly lower, as can be seen from 
Fig. 13 and Tab. 13. 
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M(VZ) 




I 1 < < I 1 < < I < < l' ' , 1 1 I , 1 1 ii I n— II 

600 800 1000 1200 1400 1600 1800 

M(VZ) [GeV] 



Figure 12. Distributions of the mass of the VZ-system for the £^£ + Aj channel with full set of 
cuts, Tab. 3 and Tab. 18. 



Mcut 


no-Higgs 


KM-LET 


JAM A 


JAM B 


lAM D 


SM 


Z + Aj 


600 


.308 


.274 


.532 


.388 


.795 


.0969 


.269 


800 


.204 


.173 


.369 


.294 


.424 


.0500 


.128 


1000 


.120 


.0975 


.138 


.223 


.211 


.0258 


.059 


1200 


.0657 


.0511 


.0396 


.167 


.156 


.0154 


.027 



Table 12. Total cross section (in fb) for the £'^£ + 4j channel with the full set of cuts in Tab. 3 
and Tab. 18 in function of the minimum jcjc^£ invariant mass (in GeV). 



PDFs 



PDFs 





Figure 13. PDFs of the number of events in the £+£-+Aj channel for M„tt = 600 GeV, L = 50 fb"^ 
(left) and A4„t = 800 GeV, L = 200 fb"^ (right). 



The pattern of the distributions in Fig. 13 is analogous to the one in the 4:jlv channel: 
models with light resonances are easily distinguished from the SM and models without 
resonant states give smaller predictions which are in rough agreement among themselves. 
With L = 200 fb~^ all models have a high probability of producing a number of events 
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Li fb-i) 


Mcut (GeV) 


no-Higgs 


KM-LET 


lAM A 


lAM B 


lAM D 


50 


600 


68.36% 


57.90% 


97.43% 


85.61% 


99.95% 


200 


800 


97.70% 


92.81% 


99.99% 


99.97 .% 


100% 



Table 13. PBSM@95%CL, defined in Sub-Sect. 4.1, for the £+£- + 4j channel for L = 50 fb"^ 
and L = 200 fb"^ 

larger than the 95%CL for the SM. The KM-LET model, which is the one least likely to 
be detected sports a PBSM@95%CL of 93%. At L = 50 fb"^ the PBSM@95%CL for the 
KM-LET case drops to about 60%. 

4.3.3 The 3£iy + 2j channel 

This final state does not suffer from the huge backgrounds stemming from ti or V + Aj 
production. The final state corresponds to ZW production in the underlying scattering 
process. Hence scalar resonances cannot be formed. This is clearly visible in Fig. 14 where 
we show on the left the boson boson invariant mass distribution and on the right the PDF 
for L = 200 fb~^ and Mcut = 600 GeV. One notices immediately, comparing with the 
results obtained for the two previous channels, that the PDF for the lAM A model, which 
develops a scalar resonance, is now closer to the SM predictions than the KM-LET model. 




M(ZW) evts 



Figure 14. Distributions of the mass of the W Z-sysiem for the Mv + 2j channel with full set of 
cuts, Tab. 3 and Tab. 18 (left). PDFs of the number of events for Mc„t = 600 GeV and L = 200 fb~^ 
(right). 

The total cross sections as a function of the boson boson invariant mass lower cut are 
given in Tab. 14. The reported values are now in attobarns and this is reflected in a small 
expected number of events at the LHC. The PBSM@95%CL are shown in Tab. 15. With 
50 fb~^ of data only the I AM B and lAM D models have a reasonable chance of manifesting 
themselves in sizable deviations from the SM. At L = 200 fb~^ these two models would be 
almost certainly revealed while the PBSM@95%CL for both the no-Higgs and the KM-LET 
models are above 80%. 



-26- 



Mcut 


no-Higgs 


KM-LET 


JAM A 


lAM B 


lAM D 


SM 


600 


53.8 


51.2 


33.5 


101. 


92.6 


17.0 


800 


37.6 


34.3 


19.2 


84.6 


78.8 


8.92 


1000 


25.6 


22.5 


11.0 


72.0 


68.7 


4.42 


1200 


17.3 


14.7 


6.21 


57.1 


56.1 


2.38 



Table 14. Total cross section (in ab) for the Siv + 2j channel with the full set of cuts in Tab. 3 
and Tab. 18 in function of the minimum Siv invariant mass (in GeV). 



L{ fb-i) 


M,ut (GeV) 


no-Higgs 


KM-LET 


lAM A 


lAM B 


lAM D 


50 


600 


46.30% 


43.27% 


21.51% 


83.71% 


79.51% 


200 


600 


84.12% 


80.81% 


43.47% 


99.73% 


99.38% 



Table 15. PBSM@95%CL, defined in Sub-Sect. 4.1, for the 3iiy + 2j channel. 



4.3.4 The U + 2j channel 

The present channel has some similarities with the 3^i^ + 2j channel above. In some sense it 
is complementary to it. In fact it corresponds to the formation of a ZZ pair and this implies 
that no vector resonance can be formed in this channel. In Fig. 15 the ZZ invariant mass 
distribution is shown on the left while the PDF for L = 200 fb"^ and Mcut = 500 GeV 
are on the right. Fig. 15 clearly shows that the models which produce scalar resonances 
are well separated from the SM. On the contrary, the lAM B model, in which only a 
vector resonance is present, appears even more difficult to disentangle from the light-Higgs 
scenario than the non resonant ones. 



M(ZZ) PDFs 




Figure 15. Distributions of the mass of the ZZ-system for the A£ + 2j channel with full set of 
cuts. Tab. 3 and Tab. 18 (left). PDFs for M^^t = 500 GeV and L = 200 fb"^ (right). 

In Tab. 16 the values of the cross sections in attobarns are reported. These very 
modest cross sections imply that no stringent cut can be imposed and as a consequence the 
discriminatory power of this channel is smaller than that of all other final states presented 
in this paper, as can be seen from Tab. 17. Only at L = 200 fb~^ this channel could be 
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Mcut 


no-Higgs 


KM-LET 


lAM A 


lAM B 


lAM D 


SM 


A r\r\ 

4(J0 


A A ''7 

44.7 


41.9 


72.0 


38.5 


89.6 


25.5 


500 


35.6 


32.8 


62.0 


29.5 


79.0 


18.4 


600 


28.2 


25.4 


52.9 


22.2 


67.8 


13.5 


700 


22.2 


19.4 


44.0 


16.4 


52.0 


9.64 


800 


17.8 


14.9 


34.6 


12.2 


26.0 


7.09 


900 


14.0 


11.4 


23.6 


8.97 


10.9 


5.19 



Table 16. Total cross section in ah for the 4£ + 2] channel with the full set of cuts in Tab. 3 and 
Tab. 18 in function of the minimum At invariant mass (in GeV). 

useful for the lAM A and lAM D models. 



L{ fb-i) 


Mcut (GeV) 


no-Higgs 


KM-LET 


lAM A 


lAM B 


lAM D 


200 


500 


43.86% 


36.31% 


89.60% 


27.41% 


97.36% 



Table 17. PBSM@95%CL, defined in Sub-Sect. 4.1, for the 4£ + 2j channel. 



5 Conclusions 

We have described in some detail how Unitarized Models have been implemented in the 
PHANTOM MonteCarlo. We have then compared the predictions of this new implementation 
which makes use of full six fermion final state matrix elements calculations with those 
obtained with old EVBA methods and presented in the literature. The agreement between 
the two calculations depends rather strongly on the separation between the two tag jet, a 
variable which cannot be controlled in EVBA. The shape of the distribution of the invariant 
mass of the two vector boson system is in any case quite different, in particular for smaller 
invariant masses. 

Using the new tool we have adressed the important issue of the possibility of uncovering 
new physics signals in vector boson scattering processes at LHC. For this we have considered 
some typical examples of unitarized models. These can be divided in two sets: those which 
predict the formation of resonances, which in the cases we have examined can be scalar or 
vector in nature, and those in which no resonant state is formed. 

Our analysis has concentrated on counting experiments which look for an excess of 
boson boson scattering events compared with SM expectations in the large invariant mass 
region. Our results show that these experiments can detect new physics at 14 TeV for 
all considered scenarios using a luminosity of 200 fb~^. At 50 fb^^ instead the results 
depend on the particular channel and scenario under consideration. Since models with 
no resonance, models with scalar resonances and models with vector resonances behave 
differently in different final states we stress the importance of carefully analyzing all possible 
channels and not only the ones with largest cross section or smallest background. 
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A Selection cuts 



Processes 




2j£+£-vD 


2j£^u£^u 


4j£i. 


Aj££ 


2j3£iJ 


2jAiu 


Cuts 


{W+W-) 


{ZZ) 












|r?(£±)| < 


2.0 






2.0 




2.0 




MUfJb) > 


1000 


800 




1000 


1000 


1000 


800 


\^r]{jfjb)\ > 


4.8 


4.5 


4.5 


4.8 


4.8 


4.8 




PtUc) > 








70 


60 






Prijcjc) > 










200 






priiiy) > 








200 




200 




PTmiss > 




120 




100 








PT{£+r) > 




120 






200 


200 


100 


PT{i) > 






50 










2minpT{j) < 






120 










E{3) > 


180 














max\rj{j)\ > 


2.5 




2.5 




2.8 






\v{j)\ > 


1.3 


1.9 








1.2 




\^vivj)\ > 








0.6 


1.1 


1.5 




Ar]{£j) > 


0.8 


1.3 












AR{ej) > 


1 




1.5 










AR{Zj) > 














1 


M{£j) > 


180 














M{Vj) > 










300 






\PT{h)-Mh)\ > 


220 




150 














290 














-O.G 




-O.G 










cos{S(j)zz) < 














-0.4 


AR{e+£-) < 










1.0 







Table 18. Additional selection cuts for the various channels. All masses, momenta and energies 
are expressed in GeV. 
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